从 油 机 明 肯 循环 余热 回收 系统 的 效率 预测 模型 
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(天 津 大 学 内 燃 机 燃烧 学 国家 重点 实验 室 ， 天 津 300072) 


摘 E: 有 机 朗 肯 循环 (ORC) 技术 可 以 有 效 回收 重型 车 用 柴油 机 的 排 气 热能 。 但 在 道路 运行 过 程 中 ， 其 排 气 热源 
不 稳定 ， 冷 源 的 维持 也 需要 消耗 能 量 ， 使 得 蒸发 压力 和 冷凝 压力 的 在 线 优 化 与 控制 实施 均 存 在 严峻 挑战 。 为 此 ， 建 
立 了 ORC 余 热 回 收 系统 的 效率 预测 模型 ， 以 ORC 系 统 的 热力 学 第 一 定律 模型 为 基础 ， 抓 住 蒸发 器 阶 跃 响应 的 宏观 特 
征 , 提出 了 有 效 换 热量 的 动态 修正 模型 ， 以 及 工 况 模式 切换 的 预 判 因子 , 并 采用 简化 机 理 和 数据 建 模 相 结 合 的 方法 ， 
建立 了 冷凝 器 散热 功 耗 的 经 验 模型 ， 提 出 了 基于 效率 预测 模型 的 ORC 余 热 回 收 系统 优化 控制 架构 及 策略 。 道 路 仿真 
结果 表明 了 效率 模型 和 优化 算法 的 有 效 性 : ORC 系 统 的 有 效 做 功 时 间 达 到 94%， 蒸 汽 过 热度 控制 在 5-15K 之 间 ， 冷 
凝 器 散热 功 耗 与 发 电功率 的 比值 维持 在 20% 以 下 。 
KE: 柴油 机 ; 有 机 朗 肯 循环 ， 余热 回收 ， 效率 预测 模型 
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Efficiency Prediction Model of a Diesel-Engine ORC Exhaust Heat Recovery 


System 


YANG Can, XIE Hui, ZHOU Kui, CHENG Zhang 
(State Key Laboratory of Engines, Tianjin University, Tianjin 300072, China) 


Abstract: The ORC (organic Rankine cycle) system is regarded as one of the most potential ways for 
exhaust heat recovery of heavy-duty vehicle diesel engines. But the exhaust energy fluctuates widely 
during driving cycles, also the ORC cooling will introduce extra energy consumption. Consequently, it is a 
big challenge to optimize and control the evaporation pressure and condensation pressure in the driving 
cycle. To this end, a control-oriented ORC efficiency prediction model was completed at first: basing on a 
foundation model of the first law of thermodynamics; focusing on the dominant feature of the ORC 
system step response, an effective heat absorption dynamic model and operating modes corrected factor 
were built up; and an empirical model was established to describe the energy consumption for ORC 
cooling. Afterwards, a sequence optimization algorithm was proposed to seek optimum values of the 
evaporation pressure and condensation pressure dynamically. The validation results based on a simulation 
plat were encouraging: the effective working time percentage achieves 94%, the superheat degree is 
controlled as 5K-15K, and the cooling power consumption is kept less than 20% of the ORC power 
generation. 

Key words: diesel engine; organic Rankine cycle; waste heat recovery; efficiency prediction model 
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济 性 改善 效果 中。 

但 车 辆 在 道路 运行 过 程 中 ， 其 发 动机 基本 处 于 
和 怠速、 大 负荷 、 小 负荷 等 动态 变 工 况 运 行 状态 ， 内 
燃 机 的 排 气 流量 和 温度 也 必然 呈现 出 强烈 的 瞬 变 特 
征 ， 且 变化 范围 较 大 。 前 期 研究 7 发 现 ，ORC A 
热 回收 系统 在 循环 工 况 中 呈现 出 强烈 的 动态 行为 特 
征 ， 包 括 蒸 汽 参 数 的 波动 以 及 工 况 模式 的 切换 ， 从 
而 恶化 了 其 实际 的 节 油 效果 。 法 国 石 油 研 究 院 IFP 
的 Tonalsl 等 人 研究 指出 ORC 系统 在 动态 变 工 况 条 
件 下 的 优化 控制 是 其 应 用 的 主要 瓶颈 。 

此 外 ， 热 机 不 可 能 从 单一 热源 取 热 ， 并 使 之 完 
全 转化 为 有 用 功 而 不 产生 其 它 影 响 。 ORC RAER 
收发 动机 排 气 能 量 发 电 的 同时 ， 也 通过 冷凝 器 将 低 
温 热 量 释 放 到 发 动机 冷却 系统 当中 ， 增 加 了 冷却 系 
统 的 散热 负荷 和 散热 功 耗 。Horst TA 等 中 的 研究 表 
H, ORC 系统 散热 会 增加 冷却 水 泵 的 功 耗 ， 从 而 导 
臻 其 节 油 潜力 下 降 了 1.3%-1.8%. ORC 系统 回收 的 
余热 能 越 多 ， 其 散热 需求 也 越 大 ， 冷 却 系统 的 散热 
功 耗 也 会 越 大 。 ORC 系统 发 电 和 发 动机 冷却 系统 耗 
电 两 者 之 间 存 在 着 一 定 的 盘 亏 耦合 关系 00。 

可 以 发 现 , ORC 系统 在 回收 车 用 内 燃 机 排 气 能 
量 时 ， 会 呈现 出 强烈 的 动态 特征 和 耦合 特征 ， 使 得 
ORC 系统 在 动态 变 工 况 条 件 下 的 优化 决策 和 控制 
实施 均 存 在 严峻 挑战 。 为 此 ， 本 文 以 ORC 系统 的 
动态 特征 和 耦合 特征 为 核心 ， 冒 在 建立 面向 ORC 
系统 优化 控制 的 效率 预测 模型 ， 支 持 对 蒸发 压力 和 
冷凝 压力 进行 在 线 决 策 和 控制 。 


1 问题 描述 


1.1 研究 对 象 
图 1 为 本 文 研究 对 象 的 系统 构 型 ， 采用 ORC 
技术 回收 一 台 重 型 车 用 柴油 机 增 压 涡轮 后 的 排 气 能 


量 ; 并 为 ORC 系统 匹配 了 一 套 可 控 的 冷却 系统 。 


NAE 冷却 水 泵 


图 1 柴油 机 ORC 余热 回收 系统 的 构 型 示意 图 
Fig.1 Schematic diagram of the diesel-engine ORC exhaust heat 


recovery system 
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1.2 问题 描述 

ORC 系统 优化 和 控制 研究 的 核心 任务 是 : 在 道 
路 运行 过 程 中 ， 根 据 排 气 流量 和 排 气温 度 的 瞬 变 特 
征 ， 以 系统 效率 最 优 为 目标 ， 实 时 决策 ORC 系统 
的 运行 状态 ， 并 通过 控制 手段 加 以 实现 。ORC 系统 
的 系统 效率 为 


Ex (1) 
= f (Xe | Re) 

WP, Mynn 为 系统 效率 ，E;, 为 排 气 能 量 ，Wi 为 
ORC 系统 的 循环 净 功 ，W 为 冷却 附件 的 功率 消耗 ， 
Re 表示 排 气 边界 条 件 ，X。 表示 ORC 系统 的 运行 
状态 ,包括 蔡 发 压力 、 过 热度 、 冷 凝 压力 和 过 冷 度 

想 要 在 整个 动态 运行 历程 中 ， 时 刻 保持 ORC 
系统 工作 在 最 佳 状态 ， 获 得 最 优 的 系统 效率 ， 则 需 
要 在 每 一 个 控制 步 长 ， 根 据 排 气 边界 条 件 ， 包 括 排 
气温 度 和 排 气 流量 ， 优 化 和 控制 ORC 系统 的 蒸发 
压力 、 过 热度 、 冷 凝 压力 和 过 冷 度 。 其 中 ， 根 据 近 
期 的 研究 19 可 以 发 现 ， 过 热度 和 过 冷 度 的 最 优 取 
值 比较 简单 ， 不 会 随 着 运行 工 况 的 变化 而 变化 。 然 
而 ， 鞘 发 压力 和 冷凝 压力 的 最 优化 规律 尚 不 清晰 ， 
分 而 论 之 ， 主 要 有 如 下 三 点 原因 : 


1.2.1 动态 特性 

稳 态 研究 @ 表 明 ， 最 优 的 蒸发 夺 力 与 排 气温 度 
正 相关 。 排 气温 度 和 燕 发 压力 在 数值 上 是 一 一 对 应 
的 , 但 在 动态 运行 过 程 中 , 两 者 在 时 间 上 并 不 一 致 ， 
存在 一 定 的 时 间 差 ， 因 为 发 动机 和 ORC 两 个 系统 
转化 能 量 的 时 间 尺 度 存在 较 大 差异 ， 如 表 1 所 示 。 


表 1 发 动机 与 ORC 系统 能 量 转化 特征 的 对 比 


Tablel Comparison of energy conversion features between 


engine and ORC 
发 动机 ORC 系统 
装置 工作 方式 往复 连续 
能 量 转 化 方式 TARE 间 壁 传 热 
动态 时 间 尺 度 快 : 工作 循环 慢 : 分 钟 


发 动机 往复 式 的 工作 方式 ， 以 及 爆炸 燃烧 的 能 
量 释放 形式 ， 使 得 燃油 的 化 学 能 在 一 个 工作 循环 内 
完成 到 机 械 能 的 转化 ， 其 能 量 转 化 的 动态 时 间 尺 度 
较 快 ,为 一 个 工作 循环 。 PAT, ORC 的 工作 方式 是 
连续 的 ;发 动机 瞬 态 的 排 气 能 量 经 过 间 壁 传 热 被 
ORC 工 质 吸收 ,加 热 蒸发 ， 推 动 脱 胀 机 做 功 产 生 电 
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能 ， 系 统 的 传 热 和 流动 惯性 使 得 ORC 能 量 转 化 的 
动态 时 间 尺 度 较 大 ， 需 要 几 十 秒 其 至 数 分 钟 ， 如 表 
1 所 示 。 这 种 动态 特征 使 得 燕 发 压力 和 排 气温 度 一 
一 对 应 的 稳 态 控制 律 在 动态 运行 过 程 中 无 法 实施 ， 
需要 对 蒸发 压力 进行 在 线 优 化 。 


1.2.1 工 况 特性 

Tonal’ 等 人 研究 指出 , ORC 系统 从 冷 起 动 开 始 ， 
逐渐 被 加 热 、 升 压 ， 直 至 完全 蒸发 、 过 热 ， 到 推动 
膨胀 机 对 外 输出 做 功 ， 会 经 历 一 系列 的 动态 过 程 。 
Xielil 等 人 将 其 定义 为 四 种 基本 的 工 况 模式 , 包括 起 
动 模式 、 冲 转 模 式 、 做 功 模式 和 保护 模式 。 并 指出 
工 况 模式 的 波动 会 降低 ORC 系统 的 有 效 做 功 时 间 ， 
影响 节 油 效果 的 发 挥 。 因 此 ， 在 道路 运行 过 程 
中 , ORC 系统 优化 与 控制 更 大 的 任务 是 调节 蒸发 压 
Fi, We ORC 系统 落 入 到 保护 模式 ， 绥 解 工 况 模 
式 的 波动 ， 提 高 ORC 系统 有 效 做 功 时 间 。 因 此 ， 
调节 蒸发 压力 不 单单 是 优化 ORC 系统 在 每 一 个 控 
制 步 长 上 的 系统 效率 ， 更 重要 的 是 ， 避 免 ORC R 
统 工 况 模式 的 波动 。 


蒸汽 过 热度 线 
透 平 转速 线 


z | 过 热度 安全 阅 

b | = 
2 \ 透 平 高 效 区 
t 


\ \ 转 速 低 阅 值 


党 
透 平 转速 [RPM] 


ig 
aK 
& 

# 


运行 时 间 [S] 


图 2 ORC 系统 动态 的 工作 过 程 及 模式 定义 
Fig.2 Dynamic working process and mode definition of the 
ORC system 


1.2.1 耦合 特性 

热机 不 可 能 从 单一 热源 取 热 ， 并 使 之 完全 转化 
为 有 用 功 而 不 产生 其 它 影响 。ORC 系统 在 回收 发 动 
机 排 气 能 量 发 电 的 同时 ， 也 通过 冷凝 器 将 低温 热量 
释放 到 发 动机 冷却 系统 当中 ， 增 加 了 冷却 系统 的 散 
热 负 荷 和 散热 功 耗 。Yang00 等 人 研究 表明 ，ORC 
系统 回收 的 余热 能 越 多 ， 其 散热 需求 也 越 大 ， 冷 却 
系统 的 散热 功 耗 也 会 越 大 。 ORC 系统 发 电 和 发 动机 
冷却 系统 耗 电 两 者 之 间 存 在 强烈 的 耦合 关系 ， 如 图 
3 所 示 。 而 冷凝 压力 的 取 值 对 这 种 耦合 关系 有 着 强 
烈 的 调控 作用 ， 应 随 着 发 动机 转速 和 负荷 的 增加 而 
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增加 。 换 句 话 说， 冷凝 压力 的 最 优 取 值 也 受到 排 气 
条 件 的 强烈 影响 。 


H ORC 发 电功率 外 包 络 线 
12 一 一 一 冷凝 压力 3bar 
十 0 一 冷凝 压力 4bar 
10 一 一 冷凝 压力 5bar 
十 一 冷凝 压力 6bar 
8 十 一 冷凝 压力 7bar 
| 压力 8bar Avg 小 


冷却 系统 的 功率 消耗 [kW] 


rd 


2 4 6 8 10 12, 
ORC 系 统 的 发 电功率 [kW] 


14 


图 3 冷凝 压力 对 耦合 关系 的 调控 作用 


Fig.3 Effects of the condensing pressure on the profit-loss 


relationship 


综 上 所 述 ， 在 道路 运行 过 程 中 ， 因 为 ORC 系 
统 的 动态 特性 、 工 况 特征 和 耦合 特性 ， 使 得 蒸发 压 
力 和 冷凝 压力 的 优化 和 控制 存在 严峻 挑战 ， 而 建立 
面向 ORC 系统 优化 控制 的 效率 预测 模型 是 突破 这 
一 挑战 的 关键 。 


2 效率 预测 模型 


2.1 控制 架构 

综合 考虑 ORC 系统 的 动态 特性 、 工 况 特征 和 
耦合 特性 ， 提 出 了 ORC 系统 的 优化 控制 架构 ， 如 
图 4 所 示 。 针 对 给 定 排 气 边界 条 件 Rh ， 管 理 器 和 
效率 观测 器 协同 工作 ， 优 化 得 到 最 优 的 蒸发 压力 
Pi ;和 冷凝 压力 Pica FURRY LEE Te. 5 和 透 平 转速 
Naic7y 一 起 作为 控制 器 的 控制 目标 于。 sas 在 控制 层 
面 直 接 采 用 成 熟 的 PID 控制 算法 ， 控 制 ORC 系统 
透 平 转速 、 膨 胀 比 ， 工 质 泵 流量 ， 以 及 冷却 系统 附 
件 功率 等 ,该 架构 的 关键 之 处 在 于 建立 一 个 宏观 的 、 
具有 预测 功能 的 效率 模型 。 


- 


2.2 效率 模型 

根据 图 4， 效 率 模型 的 核心 功能 是 : 在 管理 器 
的 任意 决策 时 刻 x ， 根 据 该 时 刻 的 排 气 参数 R(x) 
以 及 ORC 系统 的 状态 参数 X (xc) ， 需 要 预测 ， 如 
果 管 理 器 给 出 决策 变量 为 Xi w (K) ， 则 在 x+1 时 刻 
ORC 系统 的 工 况 模 式 、 输 出 功率 以 及 冷却 功 耗 , 将 
系统 效率 1svw C+D 反馈 给 管理 器 。 系统 效率 7 
的 定义 式 为 公式 (2). 
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管理 器 


W, 


RC 


We 


效率 模型 


模式 预测 || 功率 预测 


图 4 ORC 余热 下 
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收 系统 的 优化 控制 架构 


Fig.4 Control framework of the ORC exhaust heat recovery system 


Wrc (7)—We (7) 
system =—* E i 6 m 


= f Rew X re, opt? Yrc) 


其 中 , En 表示 发 动机 的 排 气 能 量 , Welt) 表示 ORC 
系统 未 来 一 段 时 间 内 的 发 电量 ，W(7) 表示 冷却 系 
统 所 消耗 的 功率 ，& 表示 ORC 系统 的 工 况 模式 。 


根据 ORC 系统 的 动态 特性 ， 发 动机 工作 循环 
是 其 完成 能 量 转化 的 基本 时 间 单 位 ， 如 果 以 发 动机 
的 工作 循环 为 时 间 基 准 的 话 , ORC 系统 的 动态 特性 
非常 显著 ， 其 能 量 转 换 的 时 间 尺度 远大 于 一 个 发 动 
机 工作 循环 。 因 此 ORC 系统 的 发 电量 是 一 个 动态 
模型 ， 用 Wc(7) 表示 ， 其 中 + 表示 时 间 。 相 应 的 ， 
ORC 系统 发 电 所 需要 的 冷却 功 耗 用 W(7) 表示 。 


除 此 之 外 , ORC 系统 的 效率 模型 还 需要 具备 一 
项 重要 的 功能 ， 就 是 预测 ORC 系统 的 工 况 模式 ， 
其 对 于 系统 效率 7,,, 而 言 是 一 个 0-1 关系 。 准确 地 
预测 ORC 系统 的 工 况 模式 ， 是 避免 其 模式 波动 ， 
提高 做 功 模 式 所 占 时 间 比 例 的 关键 。 


总 之 ，ORC 系统 的 效率 模型 描述 的 对 象 包括 
ORC 系统 及 其 控制 器 ， 需 要 能 够 动态 地 预测 ORC 
系统 的 工 况 模式 <E,， > KADR W, (7) 以 及 冷却 功 
耗 凤 (rz) ， 其 建 模 思路 如 图 $ 所 示 。 


利用 蒸发 器 的 传 热 动态 效应 ， 建 立 ORC 系统 
吸 热量 的 实时 预测 模型 ， 考 虑 管理 器 的 决策 变量 以 
及 控制 器 的 执行 能 力 ， 基 于 ORC 系统 的 热力 学 模 
型 ， 预测 ORC 系统 的 输出 功率 、 散 热量 以 及 过 热 
FE; 并 根据 冷却 功 耗 模型 和 工 况 模型 ， 预 测 冷却 功 
耗 和 工 况 模式 ， 在 此 基础 上 计算 ORC 系统 的 效率 。 


在 效率 模型 的 建立 过 程 中 ， 首 先 作 如 下 两 点 假设 : 
1) ”管理 器 给 出 的 决策 变量 义 kco 被 控制 器 
完全 执行 ， 不 考虑 执行 器 的 执行 过 程 和 执 
行 能 力 。 
2) 不 考虑 膨胀 机 、 发 电机 、 工 质 泵 、 冷 却 系 
统 和 阀门 的 硬件 执行 过 程 ， 只 关注 其 输 H 
E, OREK mer REPR To na> 
TERE nker 、 冷 却 附件 功率 W 和 阀门 
开 度 Valves 。 


CO 


WOW), 
= Es. 


效率 计算 


Ds 


5 ORC 系统 效率 预测 模型 的 建 模 思路 
Fig.5 Structure diagram of the ORC system efficiency modeling 


2.2.1 有 效 吸 热量 

ORC 系统 的 有 效 吸 热量 受到 蒸发 器 对 流 换 热 
过 程 的 控制 ， 考 虑 到 所 建 效 率 模型 需要 在 控制 器 中 
实时 计算 ， 因 此 ， 其 模型 不 能 够 过 于 复杂 ， 需 要 在 
首 述 精度 和 计算 实时 性 之 间 进 行 折 中 处 理 ， 然 后 将 
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模型 中 不 准确 的 部 分 用 先进 的 控制 算法 进行 补偿 。 时 刻 对 应 的 管 壁 温度 T(x+1) 。 

因此 ， 莹 发 器 对 流 换 热 过 程 的 描述 采用 基于 集 总 参 在 计算 得 到 排 气 换 热 量 O， 之 后 ， 可 以 根据 蒸 

ee l Sy 发 器 管 壁 的 集 总 传 热 模 型 ， 计 算 ORC TARAX 
在 排 气 侧 ， 排 气 外 掠 叉 排 管束 换 热 过 程 的 控制 RAE Ope ,; ， 其 控制 方程 为 


方程 包括 ， 对 流 换 热 方 程 、 能 量 守恒 方程 和 对 数 温 


差 方程 ， 如 公式 (2)、(3) 和 (4)， dt 0 (9) 
T 
= h, A, AT (2) —s 
Or, ae ae Qre. ef = hrc, Arc, eT ep ~T rc, E) (10) 
QF. = a: (Ts, ;, 二 了 (3) 如 前 面 所 述 ， 在 控制 器 功能 执行 完美 的 假设 下 ， 液 
i 相 、 两 相 和 汽 相 的 边界 相对 稳定 ， 故 可 认为 工 质 在 
Sa ae CH 蒸发 器 中 换 热 的 平均 温度 Tuc e 受到 蒸发 压力 及 cs 
In Ex,in 
AT grou 唯一 控制 ， 蒸 发 压力 高 ， 则 平均 温度 高 ，Tncz EK 
其 中 ， ha 为 排 气 侧 的 平均 换 热 系统 ， 其 服从 茹 卡 乌 于 决策 量 PB, 的 函数 。 在 这 样 的 条 件 下 ,公式 (9) 
斯 卡 斯 (Zhukauskas) 换 热 准则 ， 在 此 基础 上 可 以 和 公式 〈10) 无 法 得 到 解析 解 ， 微 分 方程 组 联 立 求 
增加 排 气 换 热 系数 的 标定 因子 42 ; hs 表示 排 气 侧 解 的 计算 负担 是 内 燃 机 控制 器 无 法 承担 的 ， 必 须 得 
的 有 效 换 热 面积 , nl, 表示 排 气 的 质量 流量 , 7 ,和 换 一 种 角度 来 解决 这 个 问题 。 
Toy 分别 为 排 气 在 蒸发 器 入 口 处 .出口 处 与 管 壁 平 ale el » TE x +1 
O 均 温度 To 的 差 值 ，cs 为 柴油 机 排 气 的 比热容 ， 工 。 RAST RIA host), FE 
= 程 上 常 取 值 为 1.088 kJ/(kg:K)。 联 立 求解 可 以 得 到 ， Poo (+1) # Pre s(x)， 那 么 工 质 的 平均 温度 必然 会 
a _hprApx 发 生 相 应 的 变化 ， 使 得 T rc, AK +1) + Trc, E(K) 。 然而 ， 
N O ne = Cp, (T exin T ep) -e Crta ) (5) EAE ae ; Ae 
ns he ee KIY = 执 量 因为 蒸发 器 管 壁 的 蓄 热 作用 ， Flyin BE Tp 不 会 立即 
公式 (5) 可 以 解读 为 两 个 部 分 ， 最 大 换 热 量 O。 ma 发 生 改 变 。 根 据 (10) 可 知 ， 其 有 效 吸 热量 
和 排 气 侧 换 热效率 7& zx ， 分 别 可 表示 为 Onc .K+1) 会 因为 工 质 平均 温度 Thc tD NR 
a ee (6) 变 而 产生 一 个 脉冲 响应 ， 脉 冲 能 量 的 大 小 等 于 蒸发 
l 器 管 壁 因为 温度 变化 而 释放 或 者 吸收 的 能 量 。 如 果 
_hecAps Pac (K+) < Poo (0) > WAR RAS ERTL; IZ 
N px, p 5 17 ca (7) g , 
= a 之 则 吸收 能 量 ， 其 值 表 征 为 
FE 在 公式 (7) 中 ,指数 项 表征 了 排 气 的 实际 换 热 热流 
一 与 可 传递 的 最 大 热流 密度 的 比值 ， 其 可 进一步 表示 5(Q sp) = CeptepAT ep a1) 
为， 其 中 , AT, 可 以 表示 为 燕 发 压力 Pes 变化 而 导致 工 
tug le (8) 质 饱 和 温度 变化 量 AT s 的 线性 关系 ， 系 数 为 4。 © 
其 中 ，s5i 为 斯 坦 顿 数 ， 是 表征 强制 对 流 强 弱 的 一 个 。 ”通过 这 种 处 理 ， 可 将 公式 《9) 和 《10) 改造 为 
CEAP. $1 越 大 ， 发 生 于 排 气 与 蒸发 器 管 壁 之 间 Ope og = Let Oe) (12) 
ATR RAL. CRAPS RATRE i ER a 
通 流 面积 之 比 。“ 越 大 ， 排 气 与 蒸发 器 管 壁 之 间 的 HP, Ter 表示 蒸发 器 动态 传 热 的 惯性 时 间 常 数 ， 


换 热 过 程 进 行 得 越 充分 ， 换 热 效率 越 高 。 需要 说 明 表征 了 蒸发 器 的 营 热 能 力 与 工 质 换 热 能 力 的 比值 。 
ee a. 518 51 LOA A CHROME et SARA 
ee dn Tan 可 通过 激励 试验 得 到 ,并且 考虑 到 流量 和 温度 
速 变换， 因此 二 以 选 代 求解 ， 在 时刻， 可 以 采 。。 “站 江 数 和 物性 参数 的 影响 ; 必 对 激 励 直 验 得 到 的 
用 第 x 时 刻 的 管 壁 温度 T(x) 来 计算 排 气 侧 的 换 热 Eee a nan Ace 
量 O_( ,然后 ,根据 管 壁 的 能 量 守恒 重新 计算 x+1 Aroro 在 此 基础 上 ， 根 据 公式 “5) 和 公式 12) 


可 以 实时 预测 ORC 系统 蒸发 器 的 有 效 吸 热量 


Orc. gC) ° 


2.2.2 工 况 模式 因子 
TURRE, 是 一 个 0-1 变量 ， 在 有 效 吸 热量 
Orce 给 定 的 情况 下 ， 受 到 决策 变量 Pro URIZ 
制 器 调节 质量 流量 臣 。， 的 共同 影响 , 但 是 管理 器 需 
要 知道 当 决 策 变 量 Pres 被 执行 ， 且 质量 流量 驴 。， 
来 不 及 调节 时 ，ORC 系统 的 工 况 模式 5， ERA 
可 能 落 入 到 保护 模式 ， 关 注 的 是 一 种 可 能 性 。 
对 于 给 定 的 有 效 吸 热量 Cac eg ， 以 及 决策 变量 
Prc3， 计 算 当 决策 变量 Pre, a 和 控制 器 调节 质量 流 
量 驹 - ,来 不 及 影响 蒸汽 状态 时 ，ORC 系统 瞬时 的 
TWIRRE, 。 也 就 是 说 ， 在 第 < 时 刻 ，ORC 系统 
TR TER RAE Tac xc 和 质量 流量 
nko y (0) 与 第 x-1 时 刻 的 值 相同 ， 那 么 要 使 蒸发 器 
出 口 蒸 汽 的 温度 达到 燕 发 压力 Pres 对 应 的 安全 温 
度 Trc3wmn， 则 需要 的 最 小 能 量 Okc omn CO F 


Qre eff, min(K) = nc, V (x —D(hxe 3,min hre 2 (x -1)) C 13 ) 


则 工 况 模式 < 为 
EK) = if (Ore, ap O- Qc ett mnf)) (14) 


其 含义 为 ， 如 果 有 效 吸 热量 Our(o 大 于 
Qkc ep minK) ， 则 ORC 系统 工作 在 做 功 模 式 ，&, 取 
值 为 1; 反之 ORC 系统 会 进入 到 保护 模式 ，6, 取 
值 为 0。 


2.2.3 冷却 消耗 功率 

ORC 系统 发 电功率 Wxc(7) 表示 当 决 策 变量 
Prea 和 prc 在 一 个 决策 周期 内 的 被 控制 器 完全 执 
行 之 后 的 预测 值 。ORC 系统 中 工 质 泵 入 口 为 状态 1 
点 ， 如 图 6 M7R, HEMEN X Pras KERN 
Herca: 根据 工 质 泵 的 效率 wi ps 可 以 计算 得 到 2 的 
状态 ， 用 一 对 变量 〈Pucs Tree) PAS, FEHR 
为 有 ac; 3 点 的 状态 由 决策 变量 Parcs 控制 ， 用 一 
对 变量 (Proso Tres) REN Hr: 最 后 ， 
根据 膨胀 机 的 等 粒 效 率 Mec rs 可 以 计算 得 到 4 点 比 
JOA Hra TEESE LE, 调用 REFPROP 软件 ， 可 
以 计算 得 到 ORC 系统 在 该 条 件 下 的 循环 效率 为 
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Hic.a)— Aree 
H H 


H 
rci) (15) 


TRC, Cycle 


RC,3 **RC2 


则 ORC 系统 发 电功率 的 预测 值 为 
Wrc (T) = Ore, eff (7) Nec. Cycle (16) 
同时 ， 冷 凝 器 的 散热 量 为 


Cuc con (T) = Dre, off -Nre Oce) (17) 


图 6 ORC 系统 循环 示意 图 
Fig.6 Thermodynamic cycle diagram of the ORC system 


ORC 系统 冷却 消耗 功率 主要 包括 水 泵 功 耗 、 风 \ 
扇 功 耗 和 迎风 阻力 产生 的 功 耗 。 
对 于 水 泵 ， 其 功 耗 可 表示 为 
Wop = AVEAP., (18) 
AP, A, KARA EI, Ros TA 
水 泵 的 体积 流 率 ，AP,, RARTKAN LATE. SR, 
对 于 冷却 水 的 流动 循环 而 言 ， 水 泵 的 扬程 其 实 是 为 
了 克服 管 路 阻力 , 包括 沿 程 损 失 和 局 部 损失 。 因 此 ， 
水 泵 的 扬程 和 AP, 可 近似 表示 为 
Ac = 4,» (19) 
将 忽略 一 次 项 所 产生 的 误差 放 在 系数 人 中 进行 标 
定 ， 故 整合 公式 C18) 和 公式 〈19) 可 以 得 到 
Wop = 1A, bp (20) 
一 般 地 ， 冷 却 水 流量 和 散热 量 直接 相关 ， 可 以 将 公 
sk (20) 变形 为 


Wop = AAA; Oren (21 ) 
HH, Q.,, 为 冷凝 器 中 的 散热 量 。 在 公式 (21) 的 
基础 上 ， 可 以 将 4 、 帮 和 ,4 合并， 得 到 
Wop = Acp Co (22) 
同 理 ， 可 以 推导 得 到 风扇 功 耗 ， 以 及 迎风 阻力 
所 产生 功 耗 的 表征 模型 ， 如 公式 (23) 和 (24), 
Wor = Per Qua F (23) 


Wop = Yep Okaa, D (24) 


` 


RP, Opar T Opa ， 分 别 表示 风扇 冷却 和 迎风 淮 


却 从 散热 水 箱 中 带 走 的 热量 。 迎 风 冷 却 是 被 动作 用 
的 ， 受 道路 工 况 影 响 ， 而 风扇 辅助 冷却 是 可 以 主动 
调节 的 ， 但 两 者 的 效率 相差 很 大 。 如 果 是 风扇 辅助 
冷却 ， 那 么 其 功 耗 为 


Wor = BYE, AP» (25) 
如 果 是 迎风 冷却 ， 其 阻力 所 产生 的 功 耗 近似 为 
Wop = VEAP. (26) 


因为 迎风 冷却 是 直接 作用 ， 没 有 能 量 转 化 环节 ， 效 
率 近似 为 100%; 而 风扇 是 一 个 效率 较 低 部 件 ， 其 
AE RIESE 30% 以 下 。 引 入 一 个 迎风 冷却 的 比 
例 因子 w, KET ORC 系统 与 发 动机 冷却 系统 的 
集成 程度 。 当 迎风 冷却 所 占 比例 增高 时 ，ORC 放 热 
过 程 的 能 效 比 增 大 。 将 公式 (23) 和 (24) 合并 ， 
得 到 


Wo = Jo)O (27) 
根据 能 量 守 恒 ， 冷 凝 器 的 散热 量 和 散热 水 箱 的 
caren 


Qre, con = Q oon = Q raa (28) 
根据 公式 (22) 和 (27)， 可 以 得 到 放 热 过 程 

总 功 耗 的 特征 模型 
Wo = (Ae + f(@g NO nc. con (29) 
定义 公式 (29) PFU “A+ fl@,)” WEZ 
为 ORC 冷却 系统 的 能 效 系数 COP,,, ， 受 冷凝 压力 
空 制 。 根 据 YangH9 等 人 的 研究 结果 ， 当 全 部 采用 风 

扇 冷 却 时 ， 
COP, =-3675x Pic ,+68420x Peg 


Corr ( 3 0 ) 
-285200 已 ,+429300 


2.2.4 效率 模型 集成 
根据 效率 模型 的 建立 过 程 ， 将 公式 〈12) 采用 
后 向 离散 处 理 可 以 得 到 


A 
Qre, eff Ge) 一 (一 ee eff (x-1) 


(31) 
: = (Op, (£) + (Opp (E) 


CHT 
其 中 ，Ar 为 决策 步 长 ， 第 x 时 刻 排 气 的 放 热 量 
Qa O A 


Qr (£) = el By EOT exin) =T lK- D)N i E (32) 


APF, AR AC ae E BEF 2 Ui BERD 8 -1 A 
度 值 ， 而 第 < 时 刻 的 温度 值 可 由 公式 〈33) 进行 计 


$ 


CgpMgp (T Ep(K)—T (£ —1)) 
= Qg.(K)— Qkc ep (£) 
第 x 时 刻 , 因为 决策 变量 P_ (K) 变化 导致 蒸发 器 管 
壁 热 容 吸收 或 者 释放 的 脉冲 能 量 5(Ouw(x) 为 
AGAR) Stig AAT s(K) (34) 
RP, Aes) 表示 第 x 时 刻 决策 变量 P(e) 对 
应 的 工 质 的 饱和 温度 ， 与 第 x-1 时 刻 决策 变 量 
P(x 一 ]) 对 应 的 工 质 的 饱和 温度 之 差 。 
F 根 据 (14). (16) A (29) 进行 分 别 计算 第 
x 时 刻 的 工 况 模式 &,(x) REX W gle) 和 冷却 
功 耗 Wc(x) ， 其 离散 形式 为 
E n) = if Ore oy O Qr emn O) (35) 


(33) 


中 


$ 


Wrc (©) = Qkc, ep Mac, cya (K) (36) 
1 = 3 
Wo (x) = [Qe a OU -Mc oa D] ( 37 ) 
COP or, (x) 


H 此 ， 可 计算 ORC 系统 在 第 * 时 刻 的 系统 效率 为 
Wi (i) —W.(x) 
O Ep 
总 共有 四 个 系数 需要 标定 ， 分 别 为 排 气 侧 换 热 系数 
缩放 因子 Ae ,惯性 时 间 7w 的 流量 修正 因子 Arca 
和 温度 修正 因子 Arr 以 及 蒸发 器 管 壁 热 容 吸收 或 
者 释放 的 脉冲 能 量 温差 修正 因子 Ae 。 
3 效率 模型 验证 
基于 文献 四 所 示 仿 真 平台 ， 对 所 建 效率 模型 进 
行 验证 , 并 且 应 用 到 图 4 所 示 的 控制 架构 中 ,对 ORC 
系统 进行 在 线 优化 和 管理 ， 从 应 用 效果 上 进一步 验 
证 所 建立 效率 模型 的 有 效 性 。 


jami 


Ë n) (38) 


7 System 


3.1 模型 精度 验证 

在 低 转 速 小 负荷 ， 到 高 转速 大 负荷 四 个 工 况 点 
对 效率 模型 进行 标定 和 校准 , 参数 设置 如 表 2 所 示 。 
效率 模型 的 校准 结果 如 图 7 所 示 。 


表 2 参数 设置 
Table2 Parameters Setup 
IR ESC IR GRAD 一 
工 质 R123 
蒸发 压力 3 MPa 
过 热度 10 K 
冷凝 压力 0.4 MPa 


冷却 方式 AY ORES SHB) 
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发 动机 的 功率 [kW] 
图 7 效率 模型 与 仿真 模型 的 对 标 结果 
Fig.7 Calibration of the efficiency model 


虽然 无 法 获得 足够 的 试验 数据 验证 所 建立 的 


效率 模型 ， 但 是 如 果 所 建立 的 效率 模型 在 趋势 上 完 
全 能 够 与 详细 的 仿真 模型 吻合 ， 有 潜力 代 蔡 详细 的 


TT 


车 速 [km/h] 
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仿真 模型 ， 用 于 ORC 余热 回收 系统 的 在 线 优化 与 
控制 。 
3.2 模型 应 用 验证 

将 所 建立 的 效率 应 用 到 图 4 所 示 ORC 系统 的 
优化 控制 架构 中 ， 在 线 优化 决策 茜 发 压力 和 冷凝 压 
力 ， 改 善 ORC 系统 的 动态 特性 、 工 况 特性 和 耦合 
特性 。 
车 辆 满 裁 ， 且 完全 采用 风 忆 冷却， 在 HWFET 
工 况 ， 对 所 建立 效率 模型 进行 应 用 效果 的 验证 。 需 
要 说 明 的 是 : HWFET 工 况 时 间 为 765s， 时 间 是 有 
限 的 ， 因 此 ， 需 要 控制 ORC 系统 的 热 容量 ， 和 否则 
起 动 时 间 太 长 ， 反 而 无 法 验证 控制 方法 的 有 效 性 。 
因此 ， 在 该 验证 过 程 中 ， 采 用 了 较 低 的 热 容 量 ， 为 
50ks 不 锈 钢 当量 ; 在 该 条 件 下 ， 因 为 蒸发 器 管 壁 的 
鞭 热 能 力 较 差 ， 对 排 气 扰动 的 缓冲 能 力 不 足 ， 其 对 
ORC 系统 的 蒸汽 状态 影响 较 大 , 在 这 样 比较 不 利 的 
条 件 下 ， 验 证 所 提出 的 控制 架构 和 效率 模型 ， 应 该 
更 具有 普 适 性 。 验 证 结果 如 图 8 所 示 。 
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8 ORC 系统 在 HWFET 了 


[ 况 下 的 控制 效果 


Fig.8 Control effect of the ORC system during the HWFET driving cycle 


图 8, 对 ORC 余热 回收 系统 在 线 优化 与 控 


制 的 实现 过 程 阐述 如 下 : 


1) 采用 低压 起 动 策略 ， 缩 短 ORC 系统 起 动 


2) 在 起 动 过 程 中 ， 


3) ORC 系统 在 起 动 、 冲 转 之 后 ,i 


总 体 而 言 ? 


时 间 。 
冷凝 压力 为 上 限 值 
0.8MPa, 在 台 架 运行 中 该 值 会 受到 工 质 泵 
入 口 压力 的 限制 。 该 决策 使 得 冷却 系统 不 
工作 ， 一 方面 可 以 降低 冷却 水 泵 和 冷却 风 
扇 的 功 耗 ， 另 一 方面 可 以 热量 的 散失 ， 促 
ORC 系统 快速 暖 机 和 起 动 。 


入 到 做 功 
模式 。 在 该 模式 下 ， 管 理 器 与 效率 模型 协 
同 工 作 ， 完 成 三 大 功能 : 1 是 工 况 模式 的 
管理 ， 尽 可 能 避免 ORC 系统 进入 到 保护 
模式 ，2 是 动态 效率 的 优化 ， 在 一 段 时 间 
内 ， 控 制 ORC 系统 的 蒸发 压力 与 排 气温 
度 相 适 应 ， 提 高 ORC 系统 效率 ; 3 THE 
合 关系 的 平衡 ， 管 理 冷 凝 压力 ， 优 化 调节 
ORC 系统 发 电 与 冷却 系统 功 耗 之 间 的 等 
亏 关 系 。 有 具体 实现 结果 如 | 
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a) “ 当 和 车辆 突然 减速 ， 排 气 能 量 不 足 ， 蒸 
汽 温度 下 降 ， 过 热度 也 会 下 降 ， 这 时 ， 通 
过 降低 蒸发 压力 ， 维 持 蒸汽 的 过 热度 ， 避 
免 其 在 短 时 间 内 跌 破 安全 限 值 而 进入 到 保 


护 模式 。 当 然 ， 与 此 同时 ， 控 制 器 也 会 使 
得 工 质 泵 的 流量 降低 ， 与 燕 发 压力 协同 配 
合 , 在 整个 循环 工 况 中 ，ORC 系统 的 蒸汽 
过 热度 一 直 控制 在 10K 左右 的 水 平 。 

b) “ 当 车 辆 高 速 巡航 ， 或 者 加 速 时 ， 排 气 
PETE, AYRE A ROHR, E 
这 种 情况 下 ， 蒸 发 压力 逐渐 升 高 ， 保 持 
ORC 系统 的 高 效 运行 。 

c) 针对 ORC 系统 发 电 与 冷却 系统 功 耗 
之 间 僵 亏 关系 的 平衡 问题 , 当 车 辆 加 速 时 ， 
ORC 系统 的 持续 吸 热量 和 散热 量 都 较 大 ， 
冷却 系统 负 蓓 比较 大 ， 这 时 ， 通 过 提高 冷 
凝 压力 的 控制 目标 ， 增 大 冷却 系统 的 散热 
温差 和 能 效 比 ， 降 低 冷却 系统 功 耗 ， 当 车 
辆 减速 或 者 了 巡航 时 , ORC 系统 的 持续 吸 热 
量 和 散热 量 较 小 ,冷却 系统 的 功 耗 比较 小 ， 
在 这 种 情况 下 , 降低 冷凝 压力 的 控制 目标 ， 
以 提高 ORC 系统 的 转化 效率 为 主 。 
采用 所 建立 的 效率 模型 和 控制 架 


F 


Ky, ORC 系统 的 有 效 做 功 时 间 达 到 94%， 蒸 汽 过 热 


BERRI 


HÆ 5K-15K 之 间 ， 散 热 功 耗 与 发 电功率 的 比 


值 维持 在 20% 以 下 ， 在 HWFET 工 况 中 的 节 油 效果 
FJ 2.2%—2.8%- 


4 $ 


SH, EZE ORC 系统 的 动态 特性 、 工 况 特 
性 和 耦合 特性 ， 建 立 了 ORC 系统 的 效率 预测 模型 ， 
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it 


本 文 为 突破 ORC 余热 回收 系统 在 线 优 化 控制 


并 提出 了 基于 效率 预测 模型 的 在 线 优化 控制 架构 和 


策略 ， 


仿真 结果 验证 了 效率 模型 和 控制 架构 的 有 效 


性 ， 形 成 如 下 结论 : 


1) ”柴油 机 ORC 余热 回收 系统 存在 强烈 的 动 
态 特 性 、 工 况 特 性 和 耦合 特性 ， 使 得 基于 
稳 态 标定 的 优化 控制 律 在 动态 运行 过 程 
中 无 法 实施 ， 需 要 对 蒸发 压力 和 冷凝 压力 
进行 在 线 优 化 。 
2) ”以 ORC 系统 的 热力 学 第 一 定律 模型 为 基础 ， 抓 
主 蒸发 器 阶 跃 响应 的 宏观 特征 ， 提 出 了 有 效 换 热 
的 动态 修正 模型 ， 以 及 工 况 模 式 切换 的 预 判 因 
子 ， 并 采用 简化 机 理 和 数据 建 模 相 结合 的 方法 ， 
建立 了 冷凝 器 散热 功 耗 的 经 验 模型 ， 集 成 建立 了 
ORC 余热 回收 系统 的 效率 预测 模型 。 
3) “所 建 效率 预测 模型 在 仿真 平台 上 的 对 标 
误差 小 于 5%， 且 在 趋势 上 完全 能 够 与 详 
细 的 仿真 模型 吻合 ， 可 用 于 ORC RAE 
收 系统 的 在 线 优 化 与 控制 。 
4) HWFET 道路 仿真 结果 从 应 用 效果 上 验证 
了 所 建 效率 预测 模型 的 有 效 性 , ORC 系统 
的 有 效 做 功 时 间 达 到 94%， 蒸 汽 过 热度 控 
制 在 SK-15K 之 间 ， 散 热 功 耗 与 发 电功率 
的 比值 维持 在 20% 以 下 。 
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